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ABSTRACT:  In  this  paper  the  energy  of  detonation  i3  calcu¬ 
lated  for  a  ''gamma  law"  gas  using  two  simple  models  for  the 
hyaroaynamic  problem.  In  the  first  a  Taylor  wave  expansion 
follows  the  detonation.  In  the  second  a  pi3ton  follows  the 
detonated  gas  to  maintain  a  constant  pressure  in  it.  The 
energy  distribution  in  the  Taylor  wave  is  found  for  the  case 
gamma  equal  to  3*  It  is  shown  that  the  average  energy 
release  by  the  chemical  reaction  per  gram  is  given  by: 
eo  ■  Pl/2(y  -  l)^o  where  p^  is  the  C-J  pressure  and  is 
the  initial  loading  density  to  produce  that  pressure.  The 
value  of  gamma,  taken  as  constant,  can  be  evaluated  from 
detonation  pressure  data.  These  approximate  results  should 
be  of  use  in  gaining  an  understanding  of  the  disposition  of 
energy  in  a  detonating  charge  of  finite  dimensions. 
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THE  ENERGY  OF  DETONATION 


I.  introduction 


A  perennial  question  concerning  explosives  is  ’’how  much 
work  can  tne  explosive  perform  upon  its  environment  when  it 
is  de  y.onated?1-  Though  the  answer  is  bound  up  in  the  theory 
of  detonations  and  rarefactions  following  detonations,  it  has 
never  been  clearly  brought  out.  This  is  probably  due  to  the 
complexity  of  the  thermodynamic  hydrodynamic  prcble**  which  is 
usually  brought  to  mind  when  the  question  is  posed,  In  this 
paper  an  attempt  is  made  to  answer  the  question  by  semi- 
quantitative  solution  of  a  typical  detonation  problem  using 
simplifying  assumptions.  Details  of  the  energy  distribution 
and  a  summation  of  tne  total  energy  released  by  a  given,  quan  - 
tity  of  explosive  are  brought  out.  Two  models  are  choseni 

a)  a  onc-dimenslonal  detonation  In  a  closed  cylinder,  and 

b)  a  one-dimensional  detonation  in  a  cylinder  in  which  a  pis¬ 
ton  does  work  or,  the  explosive  at  a  rate  given  by  the  particle, 
velocity  in  the  steady  detonation.  It  is  found  that  the 
energy  released  per  gram  can  be  easily  accounted  for  in  both 
the  nensteady  process,  a),  and  the  steady  process,  b).  This 
energy  is  just  the  energy  of  explosion  at  constant  volume  as 
will  be.  shown  by  a  comparison  with  the  Hugoniot  relation. 

The  result  is  probably  not  new  and  might  have  been  antici¬ 
pated  but  the  method  should  lead  to  a  better  understanding  of 
the  energy  of  detonation.  It  is  believed  that  a  better  under¬ 
standing  of  the  detonation  energy  can  assist  us  in  improving 
the  efficiency  of  energy  transfer,  one  of  the  major  goals  of 
ordnance  research. 


II.  Tne  Equation  of  f.'ca t e 

The  ideal  polytropi*’  gas  ossumpticn  has  made  ,.t  possible 
to  develop  a  considerable  amount  of  knowledge  concerning  real 
gases  at  moderate  pressures  even  though  it  is  kno*n  that  real 
gases  go  not  nave  constant  specific  neats,  for  a  similar 
reason,  the  use  of  the  poly tropic  gas  equation  for  the  very 
non-ideal  gases  resulting  from  the  detonation  of  solid  explo¬ 
sives  is  very  attractive,  Couslr.erabJ e  background  re."  eeu 
built  ur  to  giv\.  us  confidence  ‘  r,  the  of  this  approxi¬ 
mation  which  may  be  expressed  by  the  i sen  tropic  cniai-on: 
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In  equation  1 ) ,  "A"  Is  a  function  ol  entropy.  <f  Is  a  eon- 
.  vant  which  may  be  as  large  as  3  or  more.  This  equation  has 
t:ie  virtue  of  simplifying  both  thermodynamic  and  hydrodynamic 
calculations.  Hydrodynamic  calcula cions ,  particularly  for 
the  case  "S'  -  3,  have  shown  satisfactory  agreement  with  exper¬ 
iment  over  wide  ranges  of  pressure.  The  fit  to  isentroplc 
expansion  curves  derived  from  more  detailed  calculation  with 
equations  of  state  believed  to  be  more  acceptable  Is  quite 
good,  as  can  be  seen  ip  Figure  1.  Here,  the  pressure  computed 
by  Fivhctt  and  Cowan for  a  modified  Ki3tlakowsky--Wilson 
equation  of  state  for  65/35  cyclotoT  13  seen  to  vary  almost 
as  the  cube  of  the, density.  The  measured  pressure  of. deto¬ 
nation,  e.g.,  Deal*^',  Duff  and  Houston *3),  Mallory^’'  permit 
us  to  compute  "V  and  A  at  the  C-J  state  when  equati.  •'  ±{  is 
used.  An  extension  of  equation  l)  which  becomes  a  more- 
general  equation  of  state  will  oe  adopted,  namely: 

e  -  p/(2T  -  Dfi  2) 

where  e  is  the  internal  energy  per  gram.  Equation  2)  is  true 
for  ideal  polytropic  gases  where  X  is  defined  as  the  ratio 
of  specific  heats.  Equation  l)  for  entropy,  s,  constant  fol¬ 
lows  from  equation  2)  because  of  the  thermodynamic  relation 

rte  «  T  ds  -  pdv  3) 

Certainly  equation  2)  holds  to  the  extent  that  1)  is  true  for 
the  single  Isentrope  through  the  C-J  state  provided  we  set 
e  »  0  when  •  0,  a  mos'c  reasonable  choice.  In  the  discussion 
to  follow  the  energy  defined  by  equation  2)  Is  computed  only 
along  the  curve  s  ■  s(C--J)  except  where  noted. 


III,  The  Ava I Table  Work  In  a  Detonation  (Me iel  I) 


x  »  1 


Figure  2 


3 


NAVOKD  Report  4366 


Let  us  consider  a  model  in  VfnLch  explosive  is  loaded  into 
a  rigid  (non-expanding)  cylinder  of  unit  cross-section  and  of 
unit  length.  The  density  of  loading  is  /°Q .  The  initial 
p;  ir-sure  is  take'.;  as  zero.  Let  the  explosive  column  be  con¬ 
fined  at  each  end  by  xigid  closures.  This  gives  us  a  closed 
bomb  of  constant  unit  volume.  Figure  1.  Let  us  assure  that  a 
detonation  witn  negligible  reaction  zone  is  started  at  x  *  0, 
and  proceeds  to  the  right.  There  will  be  a  sudden  rise  to  the 
detonation  pressure,  as  the  wave  progresses,  followed  by  a 
spreading  rarefaction.  At  a  time  t  ■  i/D,  the  detonation 
front  will  reach  the  right  boundary  at  x  ■  1.  At  this  instant 
all  of  the  explosive  is  detonated.  Because  the  confining  wr.1  Is 
and  cloaures  have  nov.  moved,  the  external  work  done  by  the 
products  cf  explosior;  is  zero.  Inside  the  container  the  piv’- 
sure,  density  ar-d  particle  velocity  will  be  distriouted,  as 

jhown  by  Taylor^)  in  a  simple  (or  progressive)  rarefaction 
wave.  As  shown  by  Cole(°)  the  particle  velocity  will  be  zero 
in  (he  region  0<x£  1/2.  In  the  region  l/2<x<r  the  particle 
velocity  will  be  a  linearly  function  of  distance.  Ic  will 
have  the  value  u  «  uj,  at  x  «  1,  For  the  detonation  front  the 
following  well  knowr  conditions  express  the  conservation  of 
mass,  momentum  and  energy,  respectively: 

^l/(%  *  n/(D  -  u)  4) 

n.  .  Du2  fio  5) 

el  "  eo  *  |  Pi  [V,0  -  V^l]  6> 

The  C-J  '-/potesls  imposes  the  further  .Limitation: 

D  "  U1  *  C1  7) 

is  the  sound  speed  limned  lately  behind  the  detonation  front. 
The  polytrcpic  gas  assumption  (Eq,  2)  xt'ads  to  the  result: 

D/u1  «  y  +  1  8) 

This  follows  from  equations  l),  5)  ana  7)  ana  the  defi¬ 

nition  of  sound  speeds  c^  =  (dp/d  <"*).-> 

From  4)  and  5) 

P-^  “  (D  -  U-j  )  ^ 
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With  equal  Lon 

7) 

th Lc  gives: 

Pi  *  ci^lul 

Q\ 
^  / 

By  definition 

of 

c  and  equation  l) 

c2  .  v  A  f°r-l  ■  Jf  p/C 

10) 

By  equating  pressures  in  state  l): 

Ci2^lA  =  C^Ui 

ii) 

Therefore : 

cl/ul  -  y 

12) 

and 

D/uj  «  2f  +  1;  D/c-^  =  (  y  +  l)/  jr 

13) 

The  Taylor  expansion  Is  characterized  by  a  centered  rarefac¬ 
tion  as  shown  in  Figure  3. 


Figure  3 
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In  tne  absence  of  the  left  boundary  restraint  the  expansion 
would  be  a  continuous  fan  of  straight  lines  terminating  at  a 
particle  velocity  consistent  with  c  *  0.  The  left  boundary 
condition..  ...  0  at  x  =  0,  causes  the  wave  to  cease  expanding 

when  u  -  0.  This  ervates  a  steady  region.  III.  The  straight 
lines  x/t  =  constant  are,  by  Riemann'a  theory,  (see,  e.g., 
Courant  &  Friedricks (?/  of  slope  equal  to  u  +  c  wnicn  remains 
constant  on  each  of  the  lines.  At  the  detonation  front 
u  +  c  »  Ut  +  *  D,  the  detonation  velocity.  In.  the  region 

of  expansion  u  is  a  monatonic  furction  of  c.  The  value  of  u 
Is  derived  from  the  Rlemann  condition 

du  *  dtf*  =  cd  <°/(°  lM 

For  Lsentroplc  expansion  based  on  equation  l) 

u  -  ux  «  2(c  -  cx)/(  r-  1)  15) 

Therefore: 

x/t  xu  +  cxu^+c-'-  2(c  -  l)  16) 

When  t  *=  1/D 

x  =  (u  +  c)/D;  u  =  0  17) 

Substitutions  and  simplification  of  the  result  leads  to  the 
equation: 

xx  Y/{r-  1)]  (c/Cl  -  i/y)  18) 

or 

c/cx  *  [(  JT-  1)  x  +  ij  /y  19) 

p  and  <°  are  related  to  c  by  equation'’  10)  and  l)  so  we  may 
write: 

p/Pj  »  (c/c.  )2  */(*  -1)  20) 

-  (c/c1)2/()r'  i;  21) 

This  gives  us  p  and  p  as  functions  of  x  by  using  equation 
19).  To  solve  fhe  problem  at  hand  we  also  need  u  as  a  func- 

f  1  An  r\  f*  V  T'hl  Q  1  cj  fr  <  WOP  no  I  rwf  f  <  1  7  \  Wtr  .< 

* '  ‘  -  *  *  •  * .  U  *  *  “  •  4  *  ^  "•  •  ’Lj  V  *-  v.  U  ,  v  J  * 

u/D  =  x  -  c/D 
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or: 

u/ c  2  -  (  ^  +  1)  X/  a"-  c/c^ 


which  simplifies  by  equation  19)  to: 

u/c1  «  (2  x  -l)/*'  22) 

The  et  -  -cv  In  the  cylinder  available  to  do  work  at  the  time 
the  deconation  icaSi.zz  *-ne  point  x  »  1  is  distributed  as 
kineti,.  and  potential  and  car,  f_  c-.rp^ased  as  foll^.-.s: 

Kinetic  Energy: 

K,E.  -  /  x  u|  (°x  dx  23). 

J'/i i2 

Potential  Energy  (l/2<xcl): 


Potential  Energy  (0<x<L/2);  the  steady  region: 


24) 


(P.E.) 


Ill 


25) 


The  sum  of  these  three  terns  express  the  total  energy  in  the 
system.  Integrals  can  be  evaluated  by  substitution  of  p, 
f*  and  u  as  functions  of  x.  The  integration  is  quite  simple 
for  ft  »  3  and  only  this  solution  will  be  given  here. 


Solution  for  ■  3i 


c,  p,  /o  and  u  become  functions  of  x  only  by  equations  19) , 
20 ) „  21)  and  22).  These  functions  are  (for  y' -  3): 


c y  =  c-^  (2x  +  l)/3 
px  -  p-,  (2x  +  l)3  /27 


(°x  *  <°i  (2x  +  1^3 

ux  c1  (2x  -  l)/3 
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When  x  =  1/2;  p2  =  8  pj/27 


The  total  energy  lr.  the  system  Is  therefore: 

(2x  -  l)2  (2x  +  l)dx 

ip  dx 


26) 


The  solution  of  the  two  integrals  are  easily  found.  The 
values  after  substitution  of  the  limits  are  11/24  and  63/8 
respectively.  The  energy  is  therefore: 

E  -  (11  +  63  «•  32)/432  -  p-j/4  27) 

In  equation  2?)  we  have  rubstituted  3p^  for  the  te^m 
i°  ic<£  by  using  equations  9)  and  12). 

If  we  recall  that  this  exercise  wa3  based  on  a  cylinder  of 
unit  cross-section  and  unit  length,  we  see  that  the  energy 
is  given  per  unit  volume  of  explosive.  Equation  27)  therefore 
says  that  the  energy  per  unit  volume  of  explosive  which  is 
available  to  do  work  13  p1#  the  detonation  pressure,  divided 
by  4.  We  shall  soon  show  that  when  T  differs  from  5>  the 
energy  is  p]/2(y  -  l).  The  energy  per  gram  is  therefore 
3  -  Pi  A  <°Q  for  Y  m  3  and  Pi/2*o(  jr  -  l)  for  any  other  value 
of  gamma.  An  Interesting  point  is  the  distribution  of  energy 
as  shown  Dy  equation  27).  We  see  that  11/108  or  about  10$  is 
kinetic  energy  which  exists  only  in  the  right  half  of  the 
cylinder.  G5/IO8  or  about  60$  is  potential  energy  in  the  right 
half  of  the  model.  The  remainder,  about  30$,  is  potential 
energy  in  the  left  half  of  the  'ylinuer.  The  distribution  of 
energy  is  shown  graphically  in  Figure  4.  The  distribution  for 
V  slightly  different  than  3  should  not  be  too  different. 


IV.  The  Available  Work  in  a  Detonation  (Model  II) 

A  simpler  approach  to  finding  the  energy  in  a  detonation 
is  now  considered.  As  before  we  use  a  cylinder  of  unit  cross- 
section  ana  ox  unit  iengt.i,  Figure  3.  A  piston  is  moveu  at 
velocity,  u  =  u-^,  to  confine  the  gaseous  products  of  detonation 
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Figure  5 


while  the  wave  propagates  to  the  right.  In  this  way  no 
expansion  of  the  gas  occurs  so  that  at  time  t  »  l/D  when 
the  detonation  front  is  at  x  *  1,  the  piston  will  be  at 
x  »  u/D  l/(  <Jf -f  l).  Pressure,  density  and  particle  veloc- 
it.v  will  be  uniformly  distributed  in  the  space  between  the 
piston  and  x  -  1.  Work  has  now  been  done  on  the  gas  by 
the  piston.  If  we  sum  up  the  potential  and  kinetic  energy 
in  tne  gas  and  subtract  out  the  work  acne,  the  residual 
will  be  energy  which  must  have  been  derived  from  the  explo¬ 
sion,  The  kinetic  energy  fcr  the  unit  volume  is  obviously 
<%uJ/2.  The  potential  energy  is  p1/(^T  -  l)<°^  per  gram  or 

Pl/V(  V  -  l)f°i  Tor  the  unit  volume  initially  used.  The  work 
done  by  the  piston  pushing  with  a  force  p-^  per  unit  area 
is  Pj  x  or  (p-,/{3r+  l)  when  the  detonation  wave  reaches  x  »  1, 
The  energy  wnich  must  have  come  from  the  chemical  decompo¬ 
sition  is  therefore: 

K  u:2/2  +  Pi  <V<  l)^i  "  Pi(2f*  1)  28) 


Substituting  Pi  »  ( +  l)/°0  u-j2  and  * *  *'/(#  +  l)  we  get 


E 


Pi 


2(tfr+  1)  TT+UTV  7_iT  r  +  1 


which  simplifies  to: 

E  a  Pl  /2(  r-  1) 

per  unit  of  initial  volume. 


29) 
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The  energy  release  per  gram  is: 

e  -  P!  /2(  ZT-  l)^o  30) 

It  is  interesting  to  note  that  the  result  expressed  in  equa¬ 
tion  30)  is  lust  the  term  e0  in  the  Hueoniot  relation, 
equation  6),  This  is  seen  as  follow?: 

From  equation  C) : 


e  «  e 
o 


?1  (l/^0  -  l/toJ/2 


3D 


writing  e1  «■  px/(  Y  -  l)(°i  and  /{ If  +  I '  we  ucf: 


y 

U  +  i)(  V-  tt 


-  i  +  y 

2  2TT+“lT 


which  simplifies  tr : 


e0  -  Pi/2(y-  1)^0  32) 

It  is  therefore  quite  clear  that  we  can  associate  the  energy 
available  to  perform  work  with  the  e0  term,  i.e.,  the  energy 
of  detonation  in  the  Hugoniot  equation,  irrespective  of 
whether  the  release  of  energy  is  by  a  detonation  process  or 
not.  A  uniform  constant  volume  explosion  in  which  a  wave  is 
absent  would,  for  example,  result  in  the  same  quantity  eD 
being  converted  from  chemical  energy  to  thermal  energy.  In 
this  case  <*1  so  that  e^  *  eQ  in  equation  6).  By  using 

the  general  sense  of  the  equation  of  state,  equation  2),  we 
find  that  the  pressure  for  a  constant  volume  explosion  is 
just  1/2  of  the  detonation  pressure.  This  result  which  is 
independent  of  o  agrees  quite  well  with  the  answer  obtained 
by  the  use  of  more  elaborate  equations  of  state. 


V,  Discussion 


The  result  of  the  foregoing  calculations  snow  that  a  deto¬ 
nation  gets  its  energy  to  sustain  itself  from  both  chemical 
reaction  and  the  isentropic  expansion  of  the  products  of 
detonation.  This  expansion  replaces  the  piston  on  which  work 
must  ue  done  to  sustain  a  steady  pressure  at  the  detonation 
front.  The  net  energy  released  try  uie  detonating  explosive 
is  the  same  as  that  i eleased  from  the  chemical  process  for  a 
constant  volume  explosion  as  should  be  expected. 
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c)  Christian  and  Snay,  ref.  (8) 
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It  Is  interesting  to  compare  the  energy  term,  eQ,  derived 
from  detonatLon  pressure  measurements  with  the  energy  of 
reaction  given  by  arbitrarily  writing  down  a  chemical  reaction 
for  the  explosive  decomposition.  Data  for  detonation  pres¬ 
sures  are  available  for  a  few  explosives.  These  are  tabulated 
in  Table  i  and  the  value  of  energy.  e0  computed  using  equa¬ 
tion  32)  i3  compared  With  Qy  the  energy  of  reaction  as  give^ 
by  Christian  and  Snay(°^  for  two  arbitrary  reaction  pi'ocesses. 
The  results  are  in  good  agreement  except  for  TNT.  It  should 
be  not^  that  Qv  Is  calculated  leaving  the  products  at  stand¬ 
ard  state,  that  is,  29 8°K.  e0  on  the  ocher  hand  assumes  a 
reference  state  of  p  -  0  which  implies  0°K.  A  term  must 
therefore  be  added  to  the  values  of  Qv  for  an  assumed  (ideal?) 
ga3  expansion  to  0°K.  This  will  increase  these  numt  -'’s  oy 
about  5-10$.  On  the  other  hand  the  assumption  of  constant  tf 
throughout  an  isentroplc  expansion,  (;.s  we  have  done  in 
deriving  e?)  also  tends  to  underestimate  the  value  of  e0.  The 
agreement  is  therefore  brought  into  line  by  compensating  errors, 
it  might  be  said. 

There  are  refinements  that  can  be  taken  to  arrive  at  more 
precise  calculations.  We  are  already  working  on  a  more  suit¬ 
able  equation  of  state  which,  though  arbitrary,  will  permit 
us  to  associate  one  low  pressure  points  in  the  pressure-density 
relation  aj.ong  the  isentrope  with  the  temperature.  The  usual 
fhemochemical  standard  state  can  then  be  established  for 
this  curve.  In  this  equation  will  be  a  variable  which 
approaches  the  ideal  gas  value  at  low  pressure.  Values  of  Y 
other  than  3  should  be  tried  in  the  present  analysis  to  see 
the  effect  of  2f  on  the  energy  distribution.  This  we  expect 
to  do. 

By  the  introduction  of  a  straightforward  method  to  inte¬ 
grate  the  energy  in  a  simple  detonation  model,  we  hope  to 
have  suggested  the  way  towards  the  computation  cf  efficiency 
of  energy  transfer  in  systems  which  may  be  closer  to  realis¬ 
tic. 


VI.  Conclusions 

1.  The  energy  and  the  energy  distribution  in  a  Taylor 
wave  has  been  computed  for  a  gamma  law  with  gamma  equal  to  3. 

2.  The  energy  which  must  come  from  the  chemical  reaction 
has  been  computed  for  2  models  wo] I  as  for  the  Hugoniot 
relation  and  the  results  agree.  The  total  energy  of  the 
explosion  products  above  a  ground  state  taken  at  p  =  0  is 
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found  to  be  equal  to  pj_/?( -  1)  <°  0  w-»ere  p-j_  is  the  detonation 
pressure  and  (° 0  is  the  initial  loading  density  to  produce 
this  pressure. 

3.  Energies  calculated  by  the  method  given  here  using 
measu' cd  de^cn?^ion  pressure  data  give  reasonable  agreement 
to  the  energy  of  explosion  ut  -  '"-'■'•"t-.  volume  calculated  from 
the  chemical  decomposition. 

4.  In  the  case  of  TNT,  the  energy  calculated  from  deto¬ 
nation  pressui c  is  too  low.  This  suggests  that  both  the 
chemical  reaction  and  the  measured  pressure  should  be  reexam¬ 
ined.  Possible  sources  of  discrepancy  r-re:  a)  an  error  in 
the  pressure  measurement;  b)  incomplete  reaction;  for  example, 
atomic  carbon  may  not  crystallize  quickly  thus  reducing  the 
energy;  c)  the  solid  carbon  may  not  transfer  its  energy  to 
the  gas  on  expansion  sc  that  the  computed  work  may  be  lower 
than  the  thermal  energy  released  by  the  reaction.  More 
refined  calculations  by  Fickett  lead  to  a  result  in  qualitative 
agreement  with  that  found  here. 

5.  The  results  obtained  in  this  paper  should  ne  of  assist¬ 
ance  in  understanding  the  detonation  process  in  charges  of 
finite  size. 

6.  Use  of  a  gamma  law  for  explosives  generalized  to  an 
energy  relation  which  becomes  a  caloric  equation  of  state  is 
a  practical  simplification  for  the  solution  of  thermo- 
hydrodynamic  problems  and  is  a  means  for  obtaining  answers 
which  reasonably  approximate  the  results  of  more  involved 
assumptions. 
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